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SEAWALL BOUNDARY CONDITION IN NUMERICAL 
MODELS OF SHORELINE EVOLUTION 


PART I: INTRODUCTION 


Overview 


1. This report provides potential users with a complete description of 
the method developed by Hanson and Kraus (1985) for implementing the seawall 
boundary condition in the shoreline change numerical model. Example runs are 
included so that users may test their programs. Computer programs written in 
FORTRAN 77 are given and explained for both explicit and implicit finite- 
difference numerical solution schemes. 

2. The governing principles for the seawall boundary condition are sum- 
marized in Part I. The physical basis of the seawall boundary condition is 
discussed in a general and descriptive way in Part II. Parts I and II provide 
background material and can be understood without knowledge of numerical model- 
ing. Technical details of the shoreline numerical model and implementation of 
the seawall boundary condition are given in Part III. Two example calculations 
and a discussion of numerical accuracy and efficiency are given in Part IV. 


The computer programs are described in Part V and listed in Appendix A. 


Purpose of Seawalls 


3. Chronic erosion is found along many portions of the coast of the 
United States and other coasts of the world. Coastal erosion is caused by di- 
verse factors. These include rise in mean sea level, increase in severity of 
incident waves, change in local magnitude and direction of incident waves (as 
produced, e.g., by a newly installed coastal structure), loss of sediment sup- 
ply from rivers and cliffs, and interruption of the local littoral drift by 
structures. If the cause of undesirable erosion in an area cannot be elimi- 
nated or corrected, then buildings, roads, and other resources will eventually 
become endangered, and some degree of shore protection must be undertaken. 
Chapter 1 of the Shore Protection Manual (SPM 1984) contains a detailed dis- 


cussion of the causes of coastal erosion and their remedial measures. 


4. The shore can be protected against erosion through the use of 
coastal structures, nonstructural procedures, such as beachfill, or a combina- 
tion of structures and nonstructural methods (SPM 1984; US Army Corps of Engi- 
neers 1981). In situations where extensive damage may occur because of storm 
waves and water intrusion, or where nonstructural procedures are not feasible, 
then seawalls, bulkheads, and coastal dikes are commonly constructed for beach 
erosion control and for preventing inundation. If the word "seawall" is used 
to describe any man-made or natural object which functions as a nonerodible 
barrier along the shoreline, the concept of "seawall" encompasses true sea- 
walls, coastal dikes, storm surge barriers, shore-connected breakwaters, bulk- 
heads, revetments, and rocky coastal cliffs. A coast may contain several such 
seawalls, and their presence must be taken into account when assessing the 
long-term (order of years) evolution of the shoreline. 

5. It is also necessary to estimate the impact of a proposed seawall in 
the design process for shore protection. Even a wide sandy beach cannot erode 
indefinitely; at some point in time the beach material will be exhausted, and 
permanent structures and resources will become exposed to wave action and in- 
undation. In such a situation, emergency protective measures will be taken, 
most likely by the construction of a revetment, bulkhead, or seawall. A nu- 
merical model of shoreline change must allow for the real world situation of 


the ultimate presence of a seawall. 
Seawalls and the Shoreline Change Model 


6. Numerical models provide a powerful means for making quantitative 
estimations of shoreline evolution. In particular, the so-called "one-line" 
numerical model, originating from the work of Pelnard-Considere (1954), has 
been widely applied in recent years. Kraus (in preparation) gives an anno- 
tated bibliography of the literature on one-line models. The term "one-line" 
typically refers to the shoreline; therefore, this model is often called the 
"shoreline" model. Despite the large number of applications of the shoreline 
model, representation of the action of a seawall in the model has received 
little attention. A seawall imposes a constraint, or boundary condition, on 
the solution (shoreline position) obtained with the model. 

7. The most obvious boundary condition imposed by a seawall is that 
the beach fronting the wall cannot move landward of it. Also, a seawall 


prevents the sediment contained behind it from entering the littoral system, 


thereby modifying the sand transport rate along the beach and possibly starv- 
ing the adjacent beach through the elimination of potential littoral material. 
In an extreme case, if the level of the beach in front of a seawall drops, 
waves will reflect from the wall instead of dissipating on the beach. Stand- 
ing waves can cause local scour that may temporarily increase transport along- 
shore or offshore, until a new, steeper equilibrium profile is achieved. The 
integrity of the seawall may be threatened when the beach elevation drops. 

8. In the literature, there has been very little discussion on repre- 
sentation of a seawall in the shoreline model or in other models. Essentially 
all of the work reported to date has been conducted by engineers associated 
with coastal engineering in Japan. More than 25 percent of Japan's 34,000-km- 
long (21,000-mile) coastline is protected by seawalls, coastal dikes, armor 
blocks, and similar structures (Ogawara 1983). 

9. In the early 1970's, Hashimoto et al. (1971) discussed the behavior 
of the longshore sand transport rate in front of a seawall armored by blocks. 
They recommended the longshore transport rate be set to zero if the shoreline 
reaches the seawall. Ozasa and Brampton (1980) treated the loss of berm in 
front of a seawall and devised prescriptions for introducing the action of a 
seawall in the shoreline numerical model. In essence, their procedure also 
consists of setting the longshore sand transport rate equal to zero at calcu- 
lation points where the berm has been removed and the shoreline has retreated 
to the seawall. Hanson and Kraus (1980) gave a procedure in the form of a 
simple shoreline adjustment, but this alone is unsatisfactory because it does 
not conserve sand volume. Tanaka and Nadaoka (1982) noted that the procedure 
of setting the transport rate to zero is not correct. They proposed two al- 
ternative methods, but unfortunately their methods appear to be arbitrary and 
incomplete. 

10. Recently, Hanson and Kraus (1985) have given an outline of a well- 
tested procedure for representing the action of a seawall in balance with the 
capability of the shoreline numerical model and in accordance with three gen- 
eral principles. The present report gives a complete description of their 
method. The physical reasoning behind the method is discussed in Part II. 

The principles upon which the method is based are: 


a. The shoreline in front of a seawall cannot recede landward of 
the seawall. 


Sand volume must be conserved. 


Io 


ec. The direction of sand transport alongshore must be preserved in 
accordance with the natural direction of the potential local 
transport. 


11. Although the above-listed principles are easy to understand, their 
implementation in a computer program is considerably involved, in particular, 
for b and c. The present report describes well-tested algorithms for imple- 


menting the seawall boundary condition in a general manner. 


Limitations of the Method 


12. The seawall constraint should be formulated on the same level of 
idealization as the shoreline model. Thus, it is not appropriate in the model 
to consider wave reflection and sea bottom scouring, and settling, flanking, 
and collapse of the seawall (for further discussion on one or more of these 
topics, see Sato, Tanaka, and Irie 1969; Silvester 1977; Toyoshima 1979; 
Walton and Sensabaugh 1979). It should be stressed that the procedure de- 
seribed here possesses the same limitations as well as the same advantages as 
the shoreline model. The seawall boundary condition is only valid to the ex- 
tent the shoreline model is valid. 

13. One of the most restrictive assumptions made in deriving the shore- 
line model is that the beach profile remains unchanged and moves seaward or 
shoreward in parallel to itself (an assumption of equilibrium of the profile). 
In nature, however, if a beach erodes to reach a vertical or nearly vertical 
seawall, due to wave reflection and scouring, the beach slope immediately in 
front of a seawall is expected to become steeper than the slope on the adjoin- 
ing beach without structures or steeper than the original beach before the 
seawall was built. 

14. The above discussion notwithstanding, examples can be found in the 
field of the growth and recovery of formerly eroded beaches fronting rough- 
faced sloping seawalls (Toyoshima 1979); nearly vertical seawalls (O'Brien 
1985); and even a vertical seawall (Berrigan 1985a,b). Because of an apparent 
lack of data in these cases, however, cause and effect have not been clearly 
distinguished. That is, it is not known with certainty whether the seawalls 
promoted growth of the beaches in front of them or if, e.g., sediment trans- 
port conditions changed to bring back the beaches with no relation to the sea- 
Wall, the seawall only initially functioning to protect the land behind it. A 


combination of the two scenarios is also possible. Toyoshima (1979) states 


that a rough-faced sloping and, ideally, permeable seawall will promote recov- 
ery by dissipating wave energy, similar to the functioning of a natural beach. 

15. Based on the results of their laboratory experiments, Hattori and 
Kawamata (1977) found that a necessary condition for the naturally occurring 
restoration of an eroded beach backed by a seawall is that a surf zone exist 
seaward of the wall. Essentially the same conclusion had been reached in an 
earlier laboratory study by Chestnutt and Schiller (1971). Clearly, results 
of simulations incorporating the seawall boundary condition in a shoreline 
model must be interpreted with caution. 

16. In order to account for an alongshore variation in beach slope, a 
mechanism to allow for cross-shore sand transport and a more complicated nu- 
merical scheme than that used in the shoreline model are required. Numerical 
models now exist which account for cross-shore transport in a schematic way. 
The "2-line" model of Bakker (1969) and Bakker et al. (1971), and the "N-line" 
model of Perlin and Dean (1978, 1983) are examples. Such models can, in prin- 
ciple, more realistically represent the beach slope in front of a seawall than 
ean the shoreline model. 

17. At present, however, these models, although more sophisticated than 
the shoreline model, have limitations for engineering use stemming from lack 
of knowledge of the physical mechanism of cross-shore sand transport. Numeri- 
cal instability and long computer run times are the main technical problems 
encountered. Relatively short calculation time is an appealing feature of the 
shoreline model. This feature, plus its demonstrated versatility for handling 
a wide range of boundary conditions, ensures the use of the shoreline model as 
an engineering tool in the foreseeable future. 

18. In summary for this section, to the extent that changes in beach 
eross section can be neglected in comparison to changes in beach planform, the 
shoreline model is a useful engineering tool for systematically investigating 
and estimating shoreline evolution over time periods of several months to sev- 
eral years. If seawalls are located along the coast, because of possible sig- 
nificant changes in beach cross section, particular caution should be exer- 
cised in interpreting model results. 

19. As progress is made, it will become desirable to incorporate the 
seawall boundary condition in models more sophisticated than the shoreline 
model. This task may prove to be difficult. Experience and familiarity with 
the implementation of the seawall boundary condition in the shoreline model 


should provide useful guidance. 


PART II: BACKGROUND FOR THE SEAWALL BOUNDARY CONDITION 
Action of a Seawall on a Beach 


20. There is remarkably little quantitative information available on 
the behavior of real beaches backed by seawalls. It has been long known that 
under certain wave conditions, a vertical seawall will accelerate erosion of 
the beach in front of it (see Russell and Inglis 1953; Sato, Tanaka, and Irie 
1969). Scour is the primary cause of this erosion. Sand is scoured from the 
sea bottom in front of a vertical seawall by the standing wave system produced 
by wave reflection at the wall. Any current, such as the longshore current, 
ean then transport the mobilized sand out of the area. If there is a contin- 
ued net loss of sand over a long period of time, the end result is that the 
beach in front of the seawall can no longer maintain the natural equilibrium 
profile and the beach slope will become steeper. Walton and Sensabaugh (1979) 
discuss this and other processes believed to enhance erosion of beaches backed 
by vertical or nearly vertical seawalls. 

21. On the laboratory scale, it has been amply demonstrated that a sea- 
wall does not always produce erosion when introduced in the active wave zone 
of a beach in equilibrium with the existing waves. A brief discussion will 
now be given of three experiments (Dorland 1940, Chestnutt and Schiller 1971, 
and Hattori and Kawamata 1977) performed using sand beaches in two-dimensional 
wave flumes. 

22. Dorland (1940) used moderately steep waves in an attempt to repro- 
duce storm conditions. He placed a vertical seawall at the shoreline of a 
beach which had been allowed to attain equilibrium under constant wave action, 
scooped out part of the bed in front of the seawall, and then continued 
applying the waves. In the two such experiments performed, the outer bar moved 
landward and the scooped out area partially filled with sand from the offshore. 
In a third series of runs. using three sets of wave conditions varying cycli- 
cally, Dorland similarly found that the scooped out beach partially recovered. 

23. Chestnutt and Schiller (1971) found that maximum erosion occurred 
if a seawall was placed on an equilibrium beach in a "critical" region lying 
from about 0.5 x, to 0.67 Xj, , Where x, is the width of the surf zone, as 
measured from the shoreline. When the seawall was moved to a position shore- 


ward of the critical region, the profile immediately seaward of the wall began 


to accrete, i.e., the previously wave-scoured region tended to be filled. 
Chestnutt and Schiller point out that the surf zone width depends, in part, 

on the wave period. Other factors being the same, the surf zone will be wider 
for longer period waves. Therefore, whether or not a seawall will tend to 
promote erosion or accretion depends on the wave conditions, which usually 
have a marked seasonal variation. 

24. Hattori and Kawamata (1977) recorded beach profile changes on a 
laboratory beach with and without a vertical seawall. For given wave con- 
ditions, the beach was allowed to attain equilibrium before introduction of 
the seawall. Incident wave steepness was varied for a fixed location of the 
seawall relative to the initial shoreline. They found the existence of a 
surf zone to be a necessary condition for recovery of an eroded seawall-backed 
beach. This result is in agreement with the findings of Chestnutt and Schil- 
ler (1971). The existence of a surf zone implies minimum wave reflection at 
the seawall. Hattori and Kawamata also found that the restoring wave con- 
ditions for a seawall-backed beach are similar to those for a natural labora- 
tory beach without a seawall. 

25. Movable bottom laboratory experiments are difficult to interpret 
because of scale effects, and longshore processes were absent in the experi- 
ments under discussion. Nevertheless, a reasonable conclusion to be drawn 
from the aforementioned work is that an eroded beach in front of a seawall 
tends to recover when the mean water level is low, the waves have mild steep- 
ness, and a sediment supply exists in the offshore. Toyoshima (1979), O'Brien 
(1985), and Berrigan (1985a,b) give examples of prototype beaches backed by 
seawalls which have become stable or have recovered to some degree. 

26. The interaction between beaches and seawalls is far from under- 
stood. A focused and intensive field monitoring effort is definitely needed 
as a first step toward achieving quantitative understanding of the influence 
of a seawall on the shoreline and beach profile. Without data, quantitative 
understanding and numerical modeling of the processes involved will be limited 


and suspect. 


Seawall at Oarai Beach, Japan 


27. The physical picture for the seawall boundary condition formulated 


by Hanson and Kraus (1985) is based on general observations of the shoreline 


10 


change at the seawall located south of Oarai Harbor, Ibaraki Prefecture, 
Japan. A location map is given in Figure 1. Shoreline change at this site 
has been extensively documented and numerically modeled (Kraus, Hanson, and 
Harikai 1985). There are two seawalls on this north-south oriented sandy 
beach facing the Pacific Ocean. The north seawall is a continuous massive 
concrete wall 2 km (1.24 miles) long and 5 m (16.4 ft) high from base to 
crown. Portions of the north seawall at Oarai-are shown in Figures 2 and 3. 
The face of the north seawall is mildly curved outward and armor blocks have 
sometimes been placed at the foot of the wall when the beach eroded. The 
south seawall is similarly constructed and 800 m (0.5 miles) long. Beach 
change at the north seawall has mainly been studied. 

28. When the shoreline reaches the seawall, the local beach slope 
becomes slightly steeper than the typical nearshore slope on this coast (which 
itself varies between approximately 1/50 and 1/70 from the beach face to the 
wave breaker line). The change in beach slope is mild and appears to be neg- 
ligible for purposes of applying the shoreline model. No drastic alteration 
in beach characteristics occurs and the beach is exposed at low tide (Fig- 
ure 3). At high tide, when the shoreline has receded to the seawall, broken 
waves slap against the face of the wall. 

29. Although the shoreline may reach the seawall at some location, it 
has been inferred on the basis of the observed and modeled long-term shoreline 
change that sand moves alongshore through the area to be deposited adjacent to 
a large groin at Oarai Harbor (Kraus, Harikai, and Kubota 1981; Mizumura 1982; 
Kraus and Harikai 1983; Kraus, Hanson, and Harikai 1985). Since alongshore 
variations in the slope of the beach in front of the seawall are small, the 
seawall does not appreciably alter the pattern of wave breaking. A surf zone 
usually exists in front of the seawall and the capacity for waves to move sand 
alongshore is retained. Sand is transported in the direction of the wave- 
induced longshore current, and the beach in front of the seawall has been ob- 


served to periodically erode and recover. 


Idealized Seawall Boundary Condition 


30. From the observations described above, Hanson and Kraus (1985) de- 
veloped the concept of the idealized functioning of a seawall for use with the 


shoreline model. They concluded that once the shoreline reaches a seawall at 
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Figure 1. Location map for the beach and seawalls at Oarai Harbor, Japan 


Figure 2. North seawall at Oarai, Japan, May 1980 


Figure 3. South end of north seawall at Oarai, Japan, May 1980 
(Seawall face in lower right portion of photograph) 


a particular location, sand cannot originate from that area. There can be a 
net gain, but no net loss, for a beach area in contact with a seawall (since 
it is assumed in the shoreline model that the beach level does not drop below 
the water line and that the beach slope does not change). However, sand can 
move alongshore through such an area, passing into and out of its boundaries, 


according to the natural direction of transport. Sand can also be deposited 
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in front of a seawall, thus allowing the beach to recover. 

31. In the shoreline model, it would be incorrect to set the transport 
rate equal to zero at a location where the shoreline makes contact with a sea- 
wall, as done in most previous treatments. Rather, the transport rate should 
be adjusted to allow calculation cells in contact with a seawall to transfer 
sand in order to conserve total sand volume and preserve the direction of its 
transport. 

32. On real beaches, sand is not always transported in the same direc- 
tion over the full length of the beach. Changes in the direction of transport 
may be produced, for example, by longshore variations in wave direction and 
wave height as caused by refraction over an irregular bottom, or by diffrac- 
tion at structures and headlands. Therefore, at one or more areas along a 
beach, it is possible that a net amount of sand is moving out of the area. 

The ways in which this can occur, and implications for shoreline change in 
the presence of a seawall, are described in the section Model Input Require- 


ments and Boundary Conditions, in Part WIA 


14 


PART III: SHORELINE MODEL AND THE SEAWALL BOUNDARY CONDITION 
Shoreline Model Review 


33. The theory of the shoreline model originated with Pelnard-Considere 
(1954). He assumed that the beach bottom, not necessarily of planar slope, 
always remains in equilibrium and, as a consequence, moves in parallel to it- 
self down to a certain depth, herein called the depth of closure. Therefore, 
one contour, or "line," is sufficient to describe changes in beach planform. 
This line is conveniently taken as the shoreline. Pelnard-Considere did not 
develop a numerical model but did give closed-form mathematical solutions for 
certain idealized cases and verified the results through laboratory experi- 
ments. Details of the numerical formulation of the model may be found in, 
e.g., Komar (1976, 1983), Le Méhauté and Soldate (1978) and Hanson and Kraus 
(1980). 

34. The purpose of the shoreline model is to simulate long-term evolu- 
tion of the shoreline or the beach planform. The governing equation for the 
shoreline position is obtained from the continuity equation for beach sediment 
(assumed to be cohesionless sand). A predictive formula for the sand trans- 
port rate is necessary to solve the governing equation. Sand transport and 
the resultant shoreline change depend on the local wind, waves, and currents, 
beach planform, boundary conditions, and constraints such as the one produced 
by a seawall. It will be assumed here that the longshore sand transport is 
produced solely by obliquely incident waves; other transport mechanisms are 
possible, such as coastal, tidal, and wind-generated currents. 

35. In the present work, it will be sufficient to use the equation for 


the shoreline position in its most basic form: 


Vo tea rr (1) 


1 
at D ox 
where 
= shoreline position, m 
= time, s 
depth of closure, m 


= volume rate of longshore sediment transport, m3/s 


M6 eo Sect St 
i 


= distance alongshore, m 


For simplicity, only longshore transport of sand is considered. It is 
straightforward to generalize Equation 1 to formally include contributions 

for cross-shore transport, as well as sediment sources and sinks. An equation 
given by Hallermeier (1979, 1983) for a limiting depth of sand motion in terms 
of the incident wave conditions has been recommended by Kraus and Harikai 


(1983) for use as the depth of closure (see also Kraus 1984). 


Model Input Requirements and Boundary Conditions 


36. In order to solve Equation 1, three kinds of information are re- 
quired: (a) the initial location of the shoreline with respect to some coor- 
dinate system (Figure 4) in which the x-axis is oriented along the trend of 
the coast and the y-axis points offshore, (b) an expression for the longshore 
sand transport rate, Q , and (c) boundary conditions for either y or Q at 
the two lateral ends of the beach. Of these, the initial position of the 


shoreline is readily obtained or assumed. 
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LATERAL BOUNDARY 
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Figure 4. Definition sketch for coordinate system, shoreline, 
seawall, and lateral boundary conditions 


37. The longshore transport rate, Q , is usually calculated from the 
"CERC" formula (SPM 1984, Chapter 4): 


Q = K'! (# Gs sini oye (2a) 
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where 
K = dimensionless empirical coefficient (of order 0.4) 
H = significant wave height, m 
Ce = wave group velocity, m/s 
) = angle of breaking waves to the shoreline, deg 


S = ratio of sand density to water density 
a' = volume of solids/total volume 
r = conversion factor from Root Mean Square (RMS) to significant wave 
height, if necessary (equals 1.416) 
The subscript b indicates quantities at wave breaking. The group velocity 


at breaking is calculated from: 


1/2 
Y 
where 


g = acceleration of gravity, m/s* 


ratio of wave height to water depth at breaking, approximately 
equal to 0.78 


~ 
" 


38. The angle ®b5 


line. It is equal to the difference between the angle the breaking waves 


is the angle of the breaking waves to the shore- 


makes with the x-axis and the angle the shoreline makes with the x-axis: 


= Th ax) 
op ae= 00 B= tan = (4) 
where 

6, = angle of breaking waves to x-axis, deg 


39. Common lateral boundary conditions are Q=0 at an impermeable 
barrier such as a long jetty or groin, and 23Q/ax = 0 ona beach that has a 
stable (fixed) shoreline position. The latter boundary condition on Q can 
also be expressed as 3dy/at = 0 (see Equation 1). 

40. In addition to lateral boundary conditions, which are necessary to 
solve any problem, it is sometimes required to constrain the solution, i.e., 


restrict movement of the shoreline position. For example, the shoreline along 


the beach backed by a seawall cannot recede behind the wall. In this report, 
the seawall constraint is referred to as a boundary condition although it is 
not a boundary condition in a true sense. 

41. Three terms will be defined to distinguish important transport sit- 
uations which can occur at a seawall. 

Minus area (Figure 5a) 

42, The expression "minus" area (minus calculation cell in the numeri- 
cal model) is applied if, at a given time, sand is transported out of both 
sides of the area. If a minus area occurs where the shoreline has eroded to a 
seawall, then the sand transport rate must be corrected in such a manner as to 
conserve sand volume and preserve direction of transport, in order to pass in- 
formation about the lateral boundary conditions. In the method described in 
this report, transport rate corrections along the beach are made in the direc- 
tion of sediment transport, i.e., in the downdrift direction. Therefore, 
minus cells are starting points for corrections. 

Plus area (Figure 5b) 

43. If sand is moving into an area from both sides at a given time, 
this condition defines a "plus" area (plus calculation cell in the numerical 
model). The terminology "plus cell" describes the reverse situation of a 
minus cell; consequently, transport rate corrections end at plus cells (or at 
lateral boundaries). 

Regular area (Figure 5c) 

44, The most common situation is for a certain quantity of sand to en- 
ter one side of an area and for a slightly different quantity of sand to leave 
the area on the opposite side. This is called a "regular" area (regular cell 
in the numerical model). Sand volume and direction of transport must be pre- 
served whether or not there is a local net gain or net loss of material. If 
the shoreline in a regular area is in contact with a seawall, no more sand can 
leave the cell than enters it. If the converse occurs, causing the nonphys- 
ical movement of the shoreline to a position landward of the seawall, the 
transport rates must be corrected in an appropriate manner to move the shore- 
line position to the seawall. 

45, If a wide beach exists in front of a seawall, it is not necessary 
to distinguish between minus areas, plus areas, and regular areas. These 
three concepts become important only when the shoreline makes contact with a 


seawall. 


c.REGULAR AREAS 


Figure 5. Conceptual diagram show- 
ing minus, plus, and regular areas 


Explicit Numerical Model 


46. Equation 1 will be discretized using a staggered grid representa- 


tion, as shown in Figure 6. For convenience, Equation 1 is reproduced here. 


ET, We 6 Gib) 


1 
at D ax 


The x-axis, which runs parallel to the trend of the shoreline, is divided into 
N calculation cells by N+ 1 cell faces (solid vertical lines in Figure 6), 


with a general cell denoted by i. On this grid, Q-points and y-points are 
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defined alternately. Q-points define calculation cell faces and y-points lie 
at the centers of cells. Subscripts denote locations of points along the 
beach. Both Q-grid points and y-grid points are separated by a constant dis- 
tance Ax alongshore; the distance between a Q-point and an adjacent y-grid 
point is Ax/2 . Lateral boundary conditions must be specified at the ends 
of the grid, e.g., at Q, and Qy,, . Alternatively, it is possible to spec- 
ify boundary conditions at yy and Yu » or impose a condition on y at one 


end of the grid and a condition on Q at the other end. 


SHORELINE 


ys 
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Figure 6. Definition sketch for finite difference discretization 


47, For simplicity, only one seawall will be considered. Its beginning 
and ending coordinates on the x-axis are denoted by YSBEG and YSEND, re- 
spectively, as shown in Figure 4. A general y-position at the seawall is de- 
noted by ys; 


48. In a standard explicit scheme, Equation 1 is discretized as 


yi = 2B (Q; 5 ©), )+y 


i+ 1 


i (5) 
where 

B = At/(2DAx) , s/m° 
At = time step, s 


AX = Space interval, m 


20 


49. For notational convenience, a prime on a quantity will denote its 
value at the next (future) time step; an unprimed quantity is evaluated at the 
present time step. Quantities at the present time step are known. In custom- 
ary notation, the next time step is denoted by a superscript n+ 1 and the 
present time step is denoted by a superscript n . The customary notation 
will be used in certain applications to follow. 

50. For the purpose of implementing a boundary condition, or con- 
straint, the explicit model is convenient since (a) only immediately neigh- 
boring values of Q; and y,; are involved, and (b) the implementation only 
involves the present shoreline position and present transport rates; no quan- 
tities at the next time step are used. 

51. If the shoreline moves landward of the position of the seawall at a 
certain grid point, thus violating the seawall constraint, the longshore sand 
transport rate must be corrected to conserve sand volume. The (nonphysical) 
erosion, or retreat, of the shoreline to a position behind a seawall, as shown 
in Figure 7, results in a nonphysical additional transport of sand out of the 
associated calculation cell. The transport rates at the cell faces must 
therefore be corrected to prevent the shoreline from moving behind the sea- 
wall. The correction must be made with consideration of the direction of 
transport at the two faces of the particular cell violating the seawall con- 
straint. Only minus cells and regular cells may require correction. The sea- 
wall constraint is never violated at a plus cell, because the shoreline always 
advances in a plus cell. 

52. The calculation procedure is described in detail next. An overview 
is as follows. First, the transport rates along the beach are calculated in 
order to determine the transport directions and to identify minus, plus, and 
regular cells. Then, as required, corrections start at either a seawall bound- 
ary or the first minus cell encountered in the search. After the starting 
cell is corrected, corrections to regular cells are made as necessary follow- 
ing the direction(s) of the longshore transport, until either a plus cell or a 


lateral boundary is reached. This procedure is repeated at each time step. 


Correction at a 
minus cell (Figure 7a) 
53. Since correction is necessary, the shoreline position Yi isles 


behind the seawall. The general principle governing transport corrections is 


that the transport rate at a downdrift cell face should be reduced to a value 
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a. Correction at a minus cell 
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b. Correction at a regular cell 


Figure 7. Conceptual diagram showing 
shoreline and transport corrections 
at minus and regular cells 


that will place the shoreline at the seawall. In a minus cell, the transport 


rates at both cell faces are directed outward; therefore, both need to be ad- 


It does not appear that the adjustments can be specified in a unique 


Hanson and Kraus (1985) calculate corrected transport rates as equal 


proportions of the original rates, as follows (with corrected quantities de- 


noted by a superscript asterisk): 
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The logic behind Equation 6 is perhaps more clearly understood by rearranging 


terms, to give, for example, 
eR fo pea (7) 


from which it is seen that, whereas Q; 


causes the shoreline to move from Ve 
to Yi , the corrected transport rate Q; moves the shoreline from a to 
ys; (as required). 

54. By substitution of Equation 6 into Equation 5, with Q; and 
Qi44 replaced by Q, and aut , respectively, it is verified that the 
corrected shoreline position is 


y, = YS; (8) 


Since the adjustment was made through use of the continuity equation, the pro- 


cedure conserved sand volume. 


Correction at a 
regular cell (Figure 7b) 
55. With corrections at the minus cell completed, adjustments continue 


for cells on both sides, following the direction of transport. The transport 
rate at an updrift face will have previously been corrected and should not be 
corrected again. Assuming for the purpose of explanation that the transport 
rate through a particular cell is in the positive x-direction, the adjusted 


* 
downdrift transport rate Oe is obtained by setting the new position yi 


1 
equal to ys; in Equation 5, to give 


ys; = y, + 2B (a, - Oita,8) (9) 


1 1 


The corrected transport rate is then 


* YS5 re VA 


iW) Shi = Seay (10) 
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As a simple check, insertion of Qi] from Equation 10 into Equation 5 gives 


the following desired result for the corrected shoreline position: 
Va = YSa (11) 
Again, this is the mathematical statement of the shoreline constraint. 


* 
56. After Qi 


grid point to determine yi 


and v are obtained, calculation moves to the next 
+1 in similar manner. Calculation proceeds from 
cell to cell in the downdrift direction along the seawall until either a plus 
cell or the end of the seawall is encountered. If a plus cell is encountered 
calculation of the shoreline position continues without necessity for correc- 
tion until another minus point or the end of the seawall is encountered. 

57. On the other side of the original minus cell, where the transport 
rate is in the negative x-direction, analogous corrections are made to trans- 
port rates as described above, i.e., to Q; . This allows determination of 
Vi The calculation then proceeds from cell to cell in the downdrift 
direction. 

58. The programs YSEXP and CORRE, discussed in Part V and listed in Ap- 
pendix A, calculate shoreline change with the explicit numerical scheme for a 


beach backed by a seawall. 


Implicit Numerical Model 


59. Compared to the straightforward development for the explicit 
scheme, as presented in the previous subsection, representation of the seawall 
constraint in an implicit numerical scheme is extraordinarily complex. In an 
implicit scheme, values of the new Q; are solved for simultaneously, over 
the whole grid, in terms of the old Q,; and other quantities. Thus, in 
checking to determine whether the seawall constraint has been violated, the 
time level halfway between the old and new time levels is involved. In the 
explicit method, transport rates of only those cells in contact with a seawall 
need to be corrected; in the implicit scheme, correction of one cell will af- 
fect all cells downdrift (whether in front of the seawall or not) and thus all 
cells downdrift require correction. Correction of all downdrift cells in- 
creases the complexity and execution time of the computation. 


60. As already discussed, the direction of sand transport must be 


oy 


preserved when correction of the transport rate is made to satisfy the seawall 
constraint. Since, in general, the transport direction can reverse along a 
beach, in an implicit scheme the transport rate must be solved for twice, 
starting independently from each of the two lateral boundaries. This doubles 
the number of calculations performed, even if no corrections are required, and 
greatly reduces the speed advantage the implicit method normally holds over 
the explicit solution method. Kraus and Harikai (1983) discuss and compare 
the relative efficiencies of the explicit and implicit numerical schemes for 
the shoreline model without inclusion of the seawall constraint. A similar 
comparison of relative efficiency, including operation of the seawall con- 
straint, is given in the examples discussed in Part IV. 

61. The finite difference equations in an implicit scheme will be de- 
rived for calculating shoreline change in the presence of a seawall. The grid 
and notation are the same as those used in the explicit scheme, described in 
the previous subsection. As the starting point, Equation 1 is rewritten to 


give equal weight to present and future values: 

2 paigiy( 12s 420) a) 
In finite difference form, Equation 12 becomes 

Ties (GH > Ohm) 2 sos | (13) 
where 


ye, = y, +B (Q; - Q,,,) (14) 
The quantity ye; can be interpreted as the shoreline position midway between 
y; and yi ; it is known since it only contains values at the present time 
step and input data. The quantity B' = At/(2D' Ax) differs from the un- 
primed version in that it contains the depth of closure at the new time step, 
Which can be calculated from the new wave conditions. 

62. It is possible to solve Equation 13 by an iterative procedure be- 
tween the yi and the Qi , aS done for example, by Le Méhauté and Soldate 
(1978). A computationally faster approach is to express the Qi in terms 


of the Yi through linearization of Equation 2. Such a linearization is 
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expected to provide an accurate approximation under typical wave conditions, 
for which the breaking wave angle is small (less than 30). The linearization 
method was introduced by Perlin and Dean (1978) for use with the CERC formula, 
Equation 2. The method was extended by Kraus and Harikai (1983) to account for 
an additional contribution arising from a systematic change in breaking wave 
height alongshore (Ozasa and Brampton 1980), as caused, e.g., by wave diffrac- 
tion. These references should be consulted for details. The final result is 


that the transport rate at the new time step can be expressed in the form 
ees ! i] i] i] 
Qpeok jplynu oct y pores ey 


where Ey and Fe are functions of the incident wave parameters. Substitu- 
tion of Equation 13 into Equation 15 gives a tridiagonal system of equations 
for the Q: 5 A tridiagonal system can be solved by an efficient standard 
algorithm, called the double-sweep algorithm. The solution is based on the 


following recurrence relation: 


q — ? q ' 
Qi EEI=Q" |e FE? (16) 
where 
Bi 
1 a 
SS Ta (ao sal.) (17) 
i i-1 
' ! ss ' ' 
ear gle wee (YSueg ni YOu) ube Seavey (18) 
ina 1s (2s be .) 
a i-1 
v = 1 
BI = BE! (19) 


63. The solution procedure, prior to making any corrections to account 
for the seawall, is as follows: 


a. Specify a boundary condition at i = 1 in terms of EE; and 
BEE 
i 


b. Solve Equations 17 and 18 for i = 2 to N, in ascending order. 
This constitutes the first sweep. 
ec. Specify a boundary condition for Ned 
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d. Solve Equation 16 for i = N to 1, in descending order. This 
step is the second sweep through the grid. 
e. Substitute the Q; into Equation 13 to obtain the new shore- 


line positions, y; 


64. The shoreline positions thus obtained at each time step must be 
compared with the position of the seawall to determine if the seawall con- 
straint was violated. If so, then the shoreline position and associated 
transport rates must be corrected. In general, when making corrections to 
satisfy the seawall constraint, it is necessary to calculate the Q. in 
ascending order, as well as descending order, so that transport corrections 
can be made in either direction. The above procedure must be repeated by 
using a recurrence relation similar to Equation 16, but which allows calcula- 
tion of Qi from the boundary condition at i= 1. This relation has the 


form: 
a ' ! ! 
Qi = PP} Qi_, + RR} (20) 


The quantities PP! and RR! depend on PP! and RR! 
it i i+1 i+1 


These quantities are defined similarly to EE i and FES in Equations 17 and 


,» respectively. 


18, and will not be written here. Expressions for these quantities and their 
solution scheme can be found in program YSIMP, discussed in Part V and listed 
in Appendix A. 

65. The time evolution of Va in the implicit scheme is shown pictori- 
ally in Figure 8a. For comparison, the analogous picture for the explicit 
scheme is given in Figure 8b. The shoreline positions y; are assumed to be 
the same in both cases. In the implicit scheme, it is seen that both the pre- 
sent values (time level n) and the future values of Q (time level n+ 1), 
entering through ay/at in Equation 12, are used to calculate the shoreline 
change from y; to yi F Since the shoreline change rates ay/at 
are constant during the time increment At , the shoreline change over a time 
step is a straight line. 

66. The shoreline position midway (in time) between y, and y; was 
previously denoted as yc, . It is seen that y-points lie on a straight line 
between two adjacent ye-points. Hence, yc-points represent possible extremes 
in shoreline position. The important implication of this is that in the im- 
plicit scheme the seawall constraint must be formulated in terms of the yc, 


and not the Viren 
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b. Explicit scheme 


Figure 8. Schematic diagram of the time evolution of a 
representative shoreline position coordinate 


67. Given shoreline position Vio position ye; can be calculated as 


(see Figure 6a and Equation 12), 
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since there is a half time step between ye; and y; and a full time step 


between ye; and ye; . In finite difference form, Equation 21 becomes: 
ye: = 23" (Qi - Qi.) + ye; (22) 


A major goal has been achieved by arriving at Equation 22, because the seawall 
constraint must be formulated in terms of yc-points. The implementation of 
the constraint is similar to that for the explicit scheme, and only an outline 
will be given. 
Correction at a minus area 

68. As in the explicit scheme, transport adjustments start at a minus 
cell and from there are performed in the direction of transport. For the mi- 
nus cell itself, the adjustment resembles that expressed by Equation 6 and 


reads as follows: 


¥ WPe a Vee , 
Q* = THe, = HD) Qi (23a) 
il i 
Ws 2 WS. 
* = tite oe ' 
ore De AGi@x Say Boe Vike (23b) 
i i 
Substitution of these corrected values into Equation 22, and using Equa- 
tion 13, verifies that the desired result has been obtained, i.e., 
Joe 
ye; = ys; (24) 


Finally, the corresponding corrected shoreline position is computed from Equa- 


tion 13 as 
Sf se 
Wa 2 (25) 


The corrected position is thus found to lie halfway between the previous ex- 
tremal position, yc; , and the seawall. 
Correction at a regular 


cell, positive transport 
69. Corrections are made by moving in the positive x-direction. Since 


the transport rate into the cell has already been adjusted in connection with 
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the previous (updrift) cell, only the transport rate out of the cell must be 
adjusted in order to satisfy Equation 20. This equation contains information 
about the upstream boundary condition. Before any adjustments are made at 
cell face i+ 1 , Equation 20 reads 


1 pS ' * 1 
es) hq a eg 


(26) 


where QF is the corrected rate made for the previous cell. This relation 
holds unless the seawall constraint was violated. If so, then QF must be 


adjusted by setting ye: equal to ys; in Equation 22, thus giving 
4 ' # _ Q# 
ys; = 2B (Q% QF) + ye; (27) 


This is easily solved for the corrected transport rate for the downdrift cell: 


or Ya = We 


eas Sh Te ) 


The procedure used to arrive at Equations 26-28 is continued in the downdrift 


direction until either a plus cell or a boundary is encountered. 


Correction at a regular 
cell, negative transport 
70. The procedure used here for making corrections downdrift, in the 


negative-x direction (on the other side of the minus cell), is very similar to 
the procedure described immediately above. The new transport rate at cell 


face i is given by Equation 16, i.e., 
oS i] * ! 
Oper ne hOt eel (29) 
Then the corrected transport is found to be 
_- YGa 


Te a (30) 


This procedure is repeated downstream until a plus cell or a boundary is 


encountered. 
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PART IV: EXAMPLE CALCULATIONS 


General Comments 


71. Two examples are presented. These hypothetical situations demon- 
strate applications of the shoreline model with an operative seawall boundary 
condition and allow checking of user implementations of the programs given in 
Appendix A. An attempt was made to give semi-plausible examples while also 
preserving clarity. This resulted in two idealized cases for which most of 
the common structures and boundary conditions could be included. The first 
example is that of an initially straight shoreline bounded on one side by a 
jetty. The beach is protected by the combination of a detached breakwater and 
a Straight seawall segment. The second example is a curved pocket beach lying 
between two headlands and protected by a curved seawall. Hanson and Kraus 
(1985) show results of several other sample calculations. 

72. In the examples, the wave field is introduced artificially; the 
breaking wave height and breaking wave angle were fabricated "by hand" to 
achieve the desired trends in shoreline movement in order to exercise the sea- 
wall constraint algorithms. The breaking wave data are set in the subprogram 
INDATA which is given in Appendix A. Values of the time and space steps and 
other parameters are entered via FORTRAN DATA statements. The names of param- 
eters and variables closely follow the notation of the main text of this re- 
port. The important exceptions are: the angle "theta," denoted as "Z," and 
the empirical coefficient "K," denoted as "K1" in the program. 

73. Both examples can be run using either the explicit or the implicit 
numerical scheme, programs YSEXP and YSIMP, respectively, in Appendix A. In 
the latter part of these programs a calculation is made to check sand volume 
conservation. It can be verified that volume is conserved to within trunca- 
tion error. 

Stability 

74. Before proceeding to the examples, the stability properties of 
the shoreline model are briefly reviewed. It can be shown (e.g., Hanson and 
Kraus 1980; Kraus and Harikai 1983) that for small breaking wave angles and 
constant wave height, Equation 1, together with Equation 2, reduces to the 
functional form of the heat equation, the governing equation derived by 


Pelnard-Considere (1954). The accuracy and stability properties of numerical 
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schemes for solving this equation are well known. Generally speaking, numer- 
ical accuracy can be improved somewhat by taking a smaller time step for a 
given space step, assuming negligible numerical truncation error. Increased 
computer execution time is the price paid for using smaller time steps. 
Therefore, one wants to balance speed of the calculation with numerical 
accuracy. 

75. Numerical accuracy should be distinguished from "physical" accu- 
racy. Numerical accuracy is a measure of how well a finite difference scheme 
reproduces the solution of a differential equation; physical accuracy is a 
measure of how well the differential equation (and the numerical solution if 
one is employed) describes the process of interest. 

76. For an explicit scheme, there is a stringent limitation (the 
Courant condition) on the size of the largest possible time step, other vari- 
ables being held constant. For small breaking wave angles, in the present 


case this condition is 
rn od (31a) 
- 2 


2K" at (#*c,) 
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where 


The quantity R. was called the "stability parameter" by Kraus and Harikai 
(1983). Equation 31a is an adequate indicator of stability in most applica- 
tions, since breaking wave angles are usually small. The stability parameter 
gives an estimate of the numerical accuracy of the solution, with accuracy 


typically increasing for decreasing values of R. . 
Example 1: Jetty and Detached Breakwater 


77. The initial condition is shown in Figure 9a. The initially 
straight 2,000-m stretch of beach is protected by a shore-parallel, detached 
breakwater and a seawall connected to a long jetty. The seawall is set back 
7 m from the initial shoreline. The jetty is assumed to be sufficiently long 
so as to act as a complete littoral barrier. The breakwater is drawn in Fig- 
ure 9 to aid visual understanding; in actuality, it lies much farther off- 


shore. The seawall and breakwater have been constructed to prevent erosion 
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of the beach adjacent to the jetty. The beach on the far left side of the 
figure is assumed to be outside the area of influence of the structures, and 
therefore its position remains fixed. 

78. Waves arrive at the site as shown in Figures 9b, c, andd. In 
these figures, the longshore distributions of the breaking wave height and 
breaking wave angle are displayed in graphic form above the related beach 
planform. The local breaking wave height and angle are mainly controlled by 
the detached breakwater. The shoreline that would result if there were no 
seawall is indicated by a dashed line. 

79. Figure 9b shows the result of waves arriving almost normal to the 
shoreline for a period of 84 hr. Convergence of waves behind the detached 
breakwater causes a bulge, or salient, to form. The wave direction then 
changes, Figure 9c, and waves arrive obliquely from the right for an elapsed 
time of 180 hr. This results in a loss of sand on the beach next to the 
jetty. The seawall prevents the shoreline from eroding farther landward imme- 
diately next to the jetty; the price paid is that more sand is removed from 
along the front of the seawall. Finally, as shown in Figure 9d, the wave di- 
rection changes again and waves arrive obliquely from the left. The wave 
shadow zone behind the detached breakwater also shifts and the potential re- 
gion for erosion moves to the middle of the seawall. Sand returns next to the 
jetty, and an eroded sector forms at the middle of the seawall. 

80. Although differences in shoreline positions with and without the 
seawall are moderate in this example, by altering the input wave conditions 
(e.g., by increasing the difference in breaking wave angle between applied 
wave conditions) a much greater disparity in resultant shorelines can be 


generated. 
Example 2: Pocket Beach 


81. The initial shoreline configuration for this example is shown in 
Figure 10a. A curved pocket beach approximately 2 km long is bounded by two 
long headlands which contain the littoral transport. A curved seawall is lo- 
cated 4 m landward of the initial shoreline. 

82. Waves first arrive obliquely from the right side of the figure for 
126 hr to produce the planform shown in Figure 10b. As a result, beach mate- 
rial moves toward the left headland. The seawall has protected the area on 


the right side of the beach, as seen by the shoreline change that would have 
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occurred without the seawall (dashed line). The incident waves then swing 

in direction and arrive obliquely from the left for 138 hr, as seen in Fig- 
ure 10c. Sand is transported past the center of the seawall to form a wide 
beach adjacent to the right headland. The beach planform in (c) is not a mir- 
ror image of (b) because, although the waves were mirror images, the initial 
shoreline conditions were different. 

83. In Figure 10c, the seawall is protecting approximately half of the 
shore, and much of the eroded sector is still located on the right side. In- 
tuition might have suggested more erosion on the leftmost side since the more 
recent waves were from the left. However, the interaction between waves and 
shoreline is nonlinear (Equation 2, the sine dependence), and the calculated 
change is different than might be expected. Finally, almost normally incident 
waves arrive to the coast for 72 hr, to give the result shown in Figure 10d. 
The beach has essentially returned to its initial planform, Figure 10a. A 
beach again exists all along the front of the seawall. 

84. In this example, the seawall protected the beach under episodes of 
oblique wave incidence, preventing excessive landward retreat of the shore- 
line. The seawall therefore worked to promote recovery of the beach (compare 
solid and dashed lines in Figure 10d). It should be cautioned that this re- 
sult is partially an artifact of the assumption of an equilibrium (constant) 
profile. In nature, the beach profile in an eroded area would probably become 
steeper than the average beach profile; it then might take a longer duration 


of the normally incident waves to cause the beach to recover. 


Comparison of Accuracy and Efficiency of the Explicit 
Scheme and the Implicit Scheme 


85. The configuration of Example 2 was used to compare the numerical 
accuracy and efficiency of the explicit and implicit numerical solution 
schemes when operating under the seawall constraint. Although the results are 
necessarily site-dependent, experience has shown the trends to be representa- 
tive and the conclusions qualitatively correct. Kraus and Harikai (1983) gave 
a similar comparison of explicit and implicit numerical schemes for shoreline 
models without the seawall constraint. 

86. The results of the comparison are shown in Table 1. The wave input 
used was that in Figure 10b and run for 120 hr. The values of key parameters 


were the same as in the previous examples: maximum wave height H,., = 3m, 
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T = 8s, DX = 50 m, and D = 6m. In the comparison, the time step, DT, was 
varied and the reference or standard case was taken to be the explicit scheme 
with DT = 6 hr. The relative accuracy with respect to the reference result, 
Ay. , Where Ay is the change in shoreline position between final and initial 


positions, is given at three locations on the left side of the beach. 


Table 1 


Stability and Accuracy of Explicit and Implicit 
Numerical Schemes with an Operative Seawall 


Stability feds (percent) 
At eee Relative Vis 
hi Ss Execution Time Isl 1s 0 eS 20 
Explicit Scheme 
1 0.08 5.30 -0.6 -3.1 0.0 
2 0.17 2.72 -0.5 -2.8 0.0 
4 0.34 1.43 -0.2 -1.2 0.0 
6 0.51 1.00 0.0 0.0 0.0 
8 0.67 unstable 
Implicit Scheme 
6 0.51 2.24 -0.5 -3.1 0.0 
12 1.01 1.19 -0.2 -2.0 0.0 
a4 2.02 0.67 -0.7 0.0 0.0 
60 5.05 0.35 13.4 2.3 14.7 
120 10.11 0.24 22.1 -7.9 23.3 


87. The results in Table 1 are qualitatively similar to those given by 
Kraus and Harikai (1983). The explicit model is computationally faster than 
the implicit model per time step; however, larger time steps can be taken with 
the implicit model while preserving reasonable numerical accuracy, allowing a 
potential overall speed advantage. For example, the implicit model with a 
time step of 24 hr and stability parameter of 2.02 is about 30 percent faster 
than the reference explicit result, yet still has acceptable numerical accu- 
racy. Engineering judgment must be exercised on a case-by-case basis to de- 
cide if a 24-hr time step will give acceptable physical accuracy. Ina simi- 
lar comparison without a seawall, Kraus and Harikai (1983) found the implicit 
model with a 6-hr time step to be comparable in accuracy and execution time to 
the reference explicit model with the same time step. As was discussed in 
Part III, the implicit model suffers a loss in efficiency when the seawall 


boundary condition is operative. 
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PART V: EXPLANATION OF COMPUTER PROGRAMS 


General Comments 


88. Here, an explanation is given of main operations performed in four 
of the five FORTRAN programs given in Appendix A. The programs are set up to 
compute the examples presented in Part IV. The final shoreline positions cal- 
culated in the examples are given in Part IV so that user implementations of 
the programs can be checked. 

89. The programs constitute the foundation of a "1-line model" and cal- 
culate shoreline change on a beach backed by a seawall by means of either the 
explicit or the implicit numerical scheme. In order to run the programs for a 
general case, wave information is needed to calculate the longshore sediment 
transport along the beach in question. Specifically, the breaking wave height 
and angle along the beach are required. The breaking wave field must be ob- 
tained from a wave calculation program such as a refraction program or from a 
combined refraction and diffraction program if large coastal structures are 
involved. It was beyond the scope of this report to include a numerical wave 
model. The breaking wave field will also be influenced by the plan shape of 
the beach (the so-called sediment-wave interaction), which changes with time. 
Numerical wave models and their relation to the shoreline change model are 
discussed by Kraus (1983). 

90. The five subprograms are called by a main program. Input wave data 
for the examples are fabricated in subroutine INDATA. The subroutine INDATA 
is elementary and will not be discussed. The longshore sand transport rate, 
computed by means of Equation 2, is calculated in subroutines YSEXP (explicit 
solution scheme) and YSIMP (implicit solution scheme). Shoreline change in 
the presence of a seawall is computed in subroutines CORRE (explicit) and 
CORRI (implicit). These latter two routines correct both the transport rate 
and shoreline position as described in Part III. 

91. Many of the algorithms are repeated in the subroutines. Comments 
are given once for each generic type of algorithm. For clarity, the programs 
are arranged to calculate for only one continuous seawall of arbitrary length 
and configuration. They can easily be generalized to handle any number of 
seawalls. 

92. In the explanations, the names of variables and line numbers refer 


to those in the indicated programs. Line numbers in parentheses refer to the 
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explicit program version. The names of most key variables in the programs are 
the same as those used in the main text of this report. They are again de- 
fined here to make the explanation more self-contained. The programs them- 
selves contain a large number of comment statements describing the operations 


performed in distinct program segments. 
Programs YSEXP and YSIMP 


93. Lines 170-190 (150-170): These statements initialize basic pa- 
rameters. YSBEG and YSEND define the beginning and end grid points of the 
seawall (Figure 4), with YSBEG < YSEND. The grid spacing is DX (in meters) 
and the time step is DT (in hours). NTIMES specifies the number of timesteps 
and IT1 and IT2 denote timesteps when the wave data are changed in the exam- 
ples. DENOM is the value of physical quantities in the denominator of Equa- 
tion 2b, evaluated for quartz sand. K1 is the empirical coefficient (K) in 
Equation 2b. The wave period is denoted by T (seconds). 

94. Lines 250-310 (240-290): Specify initial shoreline and seawall 
positions for a straight beach and seawall. 

95. Lines 370-450 (350-430): Specify initial shoreline and seawall 
positions for a curved beach and seawall. 

96. Line 570 (530): Call in wave data and renew as specified. 

97. Line 730 (680): Calculate closure depth, DCLOS, from wave 
conditions. 

98. Lines 790-830 (780-810): These lines specify boundary conditions 
for the simple cases of a fixed beach position and an impermeable long groin 
(jetty, headland). 

99. Lines 850-960 (720-750): Calculation of the longshore transport 
rate. 

100. Lines 1080-1150: In the implicit model, in order to make correc- 
tions in both directions, a reversed double sweep is necessary. The longshore 
transport rates in the arrays Q and QQ _ should be equal; a checking proce- 
dure is provided to verify this. 

101. Lines 1270 (890): After the shoreline position is calculated, 
each y; must be checked to see if it violates the seawall constraint. The 
subroutines CORRI and CORRE are called to do the check and to correct the 


shoreline positions and transport rates as necessary. 
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102. Line 1300 (970): This program segment is an error checking calcu- 
lation to verify that sand volume was conserved. It also accounts for sand 


that may have entered the system at the boundaries. 
Programs CORRE and CORRI 


103. Subroutines CORRE and CORRI are called by YSEXP and YSIMP, respec- 
tively. They recalculate the transport rate due to the possible limited vol- 
ume of sand in front of a seawall and adjust the position of the shoreline 
accordingly. 

104. Line 200 (190): A branch is made according to whether the trans- 
port rate Q; is less than, greater than, or equal to zero. A branch is nec- 
essary because the corrections must be performed in the direction of sand 
transport. 

105. Line (200): This and similar lines correspond to Equation 5. 

106. Line 260: Corresponds to Equation 20. 

107. Lines 270-310 (260-310): If the intermediate shoreline position 
YC (for the explicit scheme, position Y) is seaward of the seawall, no correc- 


tion is necessary. If not, the downstream transport rate Q; must be cor- 


+1 
rected in order to conserve sand volume. The position YC (Y) is then set to 
the corresponding position of the seawall. 

108. Lines 540-680 (540-680): Calculate as described above, but for 
the reversed transport direction. 

109. Lines 700-800 (700-800): Corrections at a minus point are com- 
puted. Sand cannot be generated in a minus cell located at a seawall. There- 
fore, the transport rates at both cell faces are corrected so that the shore- 
line will not move landward of the seawall. 

110. Lines 820-950 (820-970): This program segment operates in the same 
manner as similar segments previously described, except that here the calcula- 
tion is done in order of decreasing index since the transport is in the nega- 
tive x-direction. Calculation starts at the point to the left (lower i-values) 
of the minus cell and continues downstream until a plus cell is encountered. 

111. Lines 970-1100: After corrections are completed for grid points 
within the domain of the seawall, the same procedure must be carried out for 
the unprotected (unstructured) parts of the beach, if any. This step is nec- 
essary for the implicit scheme, since all values of Q are solved at once. 

It is not required in the explicit scheme, for which corrections are com- 


pletely determined point by point, at the present time step. 
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APPENDIX A: COMPUTER PROGRAM LISTINGS 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
4350 
460 
470 
480 
490 
500 
S10 
520 
330 
540 
550 
360 
370 
380 
590 
600 
610 
620 
630 
640 
650 
660 
670 
660 
690 
700 
710 
720 
730 
740 


C¥ Program YSEXP calculates shoreline change according to one line 


C¥ theory, taking into account the effects of a seawall. 
INTEGER YSBEG, YSEND 
REAL K1,KAP1 
DIMENSION Y(40) ,¥YS(40),Q@(41) ,2(40) ,H(40) , Y0(40) 
DATA YSBEG/26/,YSEND/40/,DX/50./,DT/6./ 
DATA DENOM/2.362 £/,NTIMES/44/,N/40/,1T1/15/,1T2/31/ 
DATA Ki/0.12/,T/8.0/,G6/9.806/, GAMMA/0.78/,RADIUS/ 12000. / 
WRITE (%,%) "X¥XXXXXEXPLICIT CALCULATIONXX¥XXX*’ 
WRITE(¥,%) *YSBEG=’,YSBEG,’ YSEND=’ , YSEND 

Cc 

C% 

C¥ Initialize arrays 


C¥ Straight shoreline 
DQ 100 I=1,N 
Q(I)=0. 
Y(I)=0. 
100 CONTINUE 
DO 105 I=1,N 
YS(I)=-7. 
105 CONTINUE 
Q(N+1)=0. 
DCLOS=0. 
GOTO 120 
cC¥ Curved shoreline 
DO 110 I=1i,N 
BET=ASIN(FLOAT (21-1) ¥DX/RADIUS) 
Y (I) =RADIUS#¥(1.-COS(BET) ) 
YO(I)=Y(1) 
110 CONTINUE 
DO 115 I=YSBEG, YSEND 
YS(I)=Y¥(1I)-4. 
115 CONTINUE 
120 CONTINUE 
Cc% 
WRITE (¥,10) (YS(I),I=1,N) 
KAP1=K1/ (16. ¥DENOM) 
C¥ @L=longshore transport rate over open boundary 
C% Q@L=0 
DO 200 IT=1,NTIMES+1 
IF(IT.EQ@.1.QOR.IT.EQ.IT1.OR.IT.EQ.1IT2) IC=1 
C¥ Subroutine INDATA computes relevant input wave data 
C%¥ at any desired time step. 
IF(IC.EQ.1) CALL INDATA(IT,1IT1,1T2,H,Z,N,DT) 
IF(IT.EQ@.NTIMES+1) THEN 
Ic=2 
IHOURS=(IT-1)*INT(DT) 


WRITE (¥, 40) 
WRITE(¥,%) °FINAL CONDITIONS (after °,IHOURS,’ 
ENDIF 
IF(IC.GE.1) THEN 
WRITE (%, 40) 
WRITE(¥,30) (Y(I),1=1,N) 
WRITE (%, 40) 
WRITE(¥,20) (@(I),I=1,N) 
ENDIF 
IF(IC.EQ@.2) GOTO 999 
IC=0 


DCLOS=2. 28%H(1)-68.5¥(H(1)/T)##2/G 
B=DT¥3600. /(2.*¥DCLOS#DX) 
B2=2.%*B 
C% 
DO 300 I=2,N 
ZBS=Z(1I) -ATAN((Y(1I)-Y(I-1))/DX) 
Q(1I)=H(1) ¥¥2¥SQRT (G/GAMMAXH (1) ) ¥KAPL¥SIN(2¥ZBS) 


A2 


hours)’ 


750 300 CONTINUE 
760 C¥ Boundary conditions: 
770 C¥ Pinned beach 


780 Q@{1)=@(2) 

790 C¥ Groin(s) 

800 C¥ Q(1)=0. 

810 Q(N+1)=0. 

820 Cx 

830 IF (YSBEG.GE.3) THEN 

840 DG 400 I=1,YSBEG-2 

850 Y(I)=Y¥ (I) -B2¥(Q(I+1)-Q(1)) 
860 400 CONTINUE 

870 ENDIF 

880 C¥ Correction of shoreline in front of seawall if necessary 
890 CALL CORRE(YSBEG, YSEND,@,B2,Y,YS) 
900 IF(YSEND.NE.N) THEN 

910 DO 500 I=YSEND+1,N 

920 Y(I)=¥(1I) -B2¥(Q(I+1)-@Q(I)) 
930 500 CONTINUE 

940 ENDIF 

950 CX 

960 CX 


970 C¥ Error calculation (DIFF: closed boundaries, AROQUT: open boundary) 
9380 CX Q@L=@L+@A(1) 

990 200 CONTINUE 

1000 999 CONTINUE 


1010 DIFF=0. 

1020 AAREA=0. 

1030 DG 600 I=1,N 

1040 DIFF=DIFF+YO(1I)-Y(1) 

1050 AAREA=AAREA+tABS (YO(I)-Y (1) ) 

1060 600 CONTINUE 

1070 ERROR=DIFF/AAREA 

1080 CX AROUT=GL¥DT¥3600. /DCLOS-DIFFXDX 

1090 Ce ERROGR=AROUT /AAREA 

1100 C¥ 

1110 C¥#* Output Fk 

1120 WRITE (*, *) 

1130 WRITE(*,%) *LOST SAND VOLUME=’,ERROR¥100,” %’ 
1140 Cx WRITE(%,%) ’ (Q@L¥DT/D-AREA) /ABSAREA*100=’ , ERROR*¥100, ’%’ 


1150 610 FORMAT (1X, ’SEAWALL POSITION’ / (1X, 10F8. 2) ) 
1160 20 FORMAT (1X, 7LONGSHORE TRANSPORT’ / (1X, 10F8.4) ) 
1170 630 FORMAT (1X, ’SHORELINE POSITION’/(1X,10F8.2)) 
1180 40 FORMAT (//) 

1190 STOP 

1200 END 
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Ck 
C% 


C¥ for the right end element. 


Cx 
20 


C% 
C% 
C% 
Ck 


SUBROUTINE CORRE (YSBEG, YSEND,@,B2,Y,YS) 
CORRE recalculates transport rates (Q) 
volume in front of a seawall and adjusts the shoreline 
- Explicit calculation scheme. 


position as necessary 


INTEGER YSBEG, YSEND 
REAL @(41),Y(40),YS(40) 


I=YSBEG 
IF(Q(I).GT.0) THEN 


Y(I-1)=Y (1-1) -B2%(Q(1I)-Q(I-1)) 


due to limited sand 


Q@ positive: Calc of shoreline Y with correction of Q@ and Y 


as necessary. 


IF (Q(I+1).GE.0) 


IF(Y(I).LT.YS(1)) 


THEN 
Y(I)=Y (1) -B2¥(Q(I+1)-Q(I)) 


DIFF=YS(1I)-Y(T) 
Q(I+1)=Q(1+1)-DIFF/B2 


Y(I)=YS(1 
ENDIF 
I=I+1 
IF(I.EQ.YSEN 
GOTO 10 
ENDIF 
K=I 
T=I+t 
IF (I.EQ. YSEND+1 


) 


D+!) 


) TH 


THEN 


GOTO 100 


EN 


Y(I-1)=¥ (1-1) -B2%(Q(1I)-Q(I-1)) 


GOTO 100 
ENDIF 
IF(I.EQ@. YSEND) 
T=I+t1 
GOTO 30 
ENDIF 
ELSE 
K=YSBEG-1 


IF (YSBEG.EQ@. 1) 
ENDIF 


THEN 


K=1 


Q@ negative: Search for a minus point. 


IF(Q(I+1).LT.0) TH 
I=I+1 
IF (I. EQ. YSEND) 
IF (Q(I+1).LE 


EN 


THEN 
0) 


THEN 


Y(I)=Y(1I) -B2#(Q(1+1)-Q(1)) 
IF(Y(I).LT.YS(1)) 


DIFF=YS(1I)-Y(T) 


THEN 


Q(I)=Q(1) +DIFF/B2 


Y(I)=Y 
ENDIF 
GOTO 30 
ENDIF 
ENDIF 
GOTO 20 
ENDIF 


S(T) 


If absent, calc Y 


Correct @ as necessary. 


Minus point: Corr of @ out of the element if shoreline moves 


behind seawall. 


YC(I) SY¥(1) -B2¥(Q(I+1)-Q(1)) 


IF(Y(I).LT.YS(1I)) 


THEN 


AY 


750 DIFF=YS(1)-Y(1) 


760 QDIFF=Q(I+1)-@(I) 

770 Q(I)=Q(1I) -DIFF/B2¥(Q(1I) /QDIFF) 

780 Q@(I+1)=Q(I+1) -DIFF/B2¥(Q(1+1)/QDIFF) 
790 Y(I) =YS(1) 

800 ENDIF 

810 Ck 


820 C¥ Calc of Y starting from element to the left of minus 
830 C¥ point or boundary. @ is negative. 


840 Ck 

850 30 DQ 40 J=I-1,K,-1 

860 Y(J)=Y (J) -B2¥(Q(J+1)-Q(J)) 
870 IF(Y¥(J).LT.YS(J).AND.J.GE.YSBEG) THEN 
880 DIFF=YS(J)-Y(J) 

890 Q(J)=@(J)+DIFF/B2 

900 Y(J)=YS(J) 

910 ENDIF 

920 40 CONTINUE 

930 T=I+1 

940 IF(1.GE.YSEND+1) GOTO 100 
950 Ce 


960 C¥ Calc of Y starting from element to the right of minus 
970 C¥ point or boundary. @ is positive. 


980 Ck 

990 GOTO 10 

1000 100 CONTINUE 

1010 RETURN = 
1020 END 


A5 


100 C¥ Program YSIMP is an implicit version of program YSEXP and 
110 C¥ calculates shoreline change according to one line theory, 
120 C¥ taking into account the effects of a seawall. 


130 INTEGER YSBEG, YSEND 

140 REAL Ki,KAP1 

150 DIMENSION 2(40),Y(40),YS(40),Y0(40) ,@(41), YCOLD(40) ,E(40) ,F (40) 
160 DIMENSION EP (40) ,FP(40),BP(40),P(41),R(41),@Q(41),H(40) 

170 DATA YSBEG/1/,YSEND/40/,DX/50./,DT/6./ 

180 DATA DENOM/2.362 /,NTIMES/56/,N/40/,1T1/22/,1T2/45/ 

190 DATA K1/0.12/,17/8.0/,G/9.806/, GAMMA/0.78/,RADIUS/12000. / 

200 WRITE (%,%) > XHXXHXHHIMPLICIT CALCULATIONX¥¥%#%%%’ 

210 WRITE(*,¥) °>YSBEG=’,YSBEG,’ YSEND=’ , YSEND 

220 CX 


230 C¥ Initialize arrays 
240 C¥ Straight shoreline 


250 DO 100 I=1,N 

260 Q(I)=0. 

270 Y(1)=0. 

280 100 CONTINUE 

290 DO 105 I=YSBEG, YSEND 
300 YS(I)S=7. 

310 105 CONTINUE 

320 Q@(N+1)=0. 

330 DOLD=0. 

340 DCLOS=0. 


ckaia) (Ge, GOTG 120 
360 C¥ Curved shoreline 


370 DG 110 I=1,N 

380 BET=ASIN(FLOAT (21-1) *DX/RADIUS) 
390 Y (1) =RADIUS¥(1.-COS(BET) ) 

400 YO(I)=Y¥(1) 

410 110 CONTINUE 

420 DO 115 I=YSBEG, YSEND 

430 VS) SYP (80) Sha 


440 115 CONTINUE 
4350 120 CONTINUE 


460 Ck 

470 WRITE(¥,10) (YS(1I),1=1,N) 

480 KAP1L1=K1/ (16. ¥DENOM) 

490 C¥ @L=longshore transport rate over open boundary 
S00 C¥ QL=0. 

SiO C#¥ C=correction term in continuity calculation 

520 CX c=1.0 

530 DO 200 IT=1,NTIMES+1 

540 IF(IT.EQ@.1.OR.1IT.EQ.IT1.OR.IT.EG@.IT2) IC=1 


S50 C¥ Subroutine INDATA computes relevant input wave data 
560 C¥ at any desired time step. 


370 IF(IC.EQ.1) CALL INDATA(IT,1IT1,1T2,H,Z,N,DT) 
380 IF(IT.EQ.NTIMES+1) THEN 

390 Ic=2 

600 THOURS=(IT-1)¥INT (DT) 

610 WRITE (¥, 40) 

620 WRITE (%¥,%) ’FINAL CONDITIONS (after ’,IHOURS,’ hours)’ 
630 ENDIF 

640 IF(IC.GE.1) THEN 

650 WRITE (¥, 40) 

660 WRITE(#,30) (Y¥(1I),1=1,N) 

670 WRITE (%, 40) 

680 WRITE(¥,20) (Q(I),1=1,N) 

690 ENDIF 

700 IF(IC.E@.2) GOTO 999 

710 IC=0 

720 DOLD=DCLOS 

730 DCLOS=2. 28%H(1)-68.5%(H(1)/T) ¥*2/G 

740 B=DT¥3600./(2.¥DCLOS#DX) 
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1300 
1310 
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1330 
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1350 
1360 
1370 
1380 
1390 
1400 


BOLD=DT¥3600. / (2. ¥DOLDXDX) 
YCOLD(1)=Y(1)+BOLD¥(@(1)-Q@(2)) 


C¥ Boundary conditions:1 
C¥ Groin causing @(1)=0. 


E(1)=0. 
F(1)=0. 


C¥ Pinned beach as @(1)=@Q(2) 


300 


E(i)=1. 
F(1)=0. 


DG 300 I=2,N 
YCOLD(I)=Y(1)+BOLD¥(Q(1I)-Q(I+1)) 
ZS=ATAN((Y(1)-Y(I-1))/DxX) 
Z2=2.%¥2Z(1) 
PWR=H (1) ¥¥2%SQRT (G/GAMMA#H(T) ) 
EP (1) =PWR¥KAPL¥2¥COS (22) %(COS(ZS) )¥#2/DX 
FP(1)=PWR¥KAPIX¥SIN(Z2) ¥(2¥(COS(ZS) )¥#2-1.) 
BP (I)=BXEP (1) 
DEN=1.+BP(1I)¥(2.-E({I-1)) 
E(I)=BP(1I)/DEN 
F(I)=(FP(1) +EP(I)¥(YCOLD(I-1)-YCOLD(I))+BP(1I)¥#F(I-1))/DEN 
CONTINUE 


C¥ Boundary condition 2: groin 


C¥ 


400 


Q(N+1)=0. 


DO 400 I=N,1,-1 
Q@(I)=E(1) ¥Q(I+1) +F (1) 
IF(IT.E@.1) THEN 
YCOLD(1I)=Y (1) +B¥(Q(I) -Q(I+1)) 
ENDIF 
CONTINUE 


C¥¥#*¥ Reversed double sweep ¥*¥*¥ 
C¥ Boundary conditions 3: groin 


C¥ 


500 


P(N+1)=0. 
R(N+1)=0. 


DO SOO I=N,2,-1 
P(I)=BP(1I)/(1.+BP(1I)¥(2.-P(I+1))) 
R(I)=(FP(1I) +EP(I)¥(YCOLD(I-1)-YCOLD(1))+BP(T)#R(I+1))/ 
(1.+BP(I)¥(2.-P(I+1))) 
CONTINUE 


C¥ Boundary condition 4 (alt 1: closed boundary, alt 2: open) 


Ck 
C% 


BA(1)=0. 
QQ(1)=R(2)/(1.-P(2)) 


DO S50 I=2,N+1 

Q@@(I)=P(1I)*¥@Q(I-1)+R(1) 

CHECK=ABS (Q@Q@(1)-Q(I)) 

IF (CHECK.GT.0.0005) WRITE(%,%) *TRANSPORT CALC. DIFFER’ 
CONTINUE 


C¥Correction of shoreline in front of seawall if necessary 


CX 
Cx 


CALL CORRI(YSBEG, YSEND,@,B,YCOLD,E,F,P,R,Y,YS,N) 


C¥ Error calculation (DIFF: closed boundaries, AROUT: open boundary) 


600 


IF(IT.EQ.NTIMES) C=0.5 

@L=QL+C¥Q(1) 

CONTINUE 

CONTINUE 

DIFF=0. 

AAREA=0. 

DO 600 I=1,N 
DIFF=DIFF+YO(1)-Y(1) 
AAREA=AAREAtABS (YO(I)-Y(1I)) 

CONTINUE 


AT 


1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 


ERRGR=DIFF/AAREA 
C¥ ARGUT=GL¥DT¥3600. /DCLOS-DIFF#DX 
C% ERRGR=AROUT /AAREA 
Cx 
C¥#¥*¥ Output Xt 
WRITE (%, ¥) 
WRITE (¥,%) "LOST SAND VOLUME=’,ERROR¥100,* %’ 
C% WRITE (¥,¥) ’(QL¥DT/DCLOS) /ABS(AREA) ¥100=’, ERROR*100, ’ 


10 FORMAT (1x, >SEAWALL POSITION’/(1X,10F8.2)) 
20 FORMAT (1X, ’LGNGSHORE TRANSPORT’/ (1X, 10F8.4)) 
30 FORMAT (1X, SHORELINE POSITION’/ (1X, 10F8.2)) 
40 FORMAT (//) 

STOP 

END 
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Cc 
C*¥ 


Ck 
C* 


C%¥ for right end element. 


C% 


20 


C% 
C¥ 
C% 
C% 


SUBROUTINE CORRI(YSBEG, YSEND,@,B, YCOLD,E,F,P,R,Y,YS,N) 
CORRI recalculates transport rates 
volume in front of a seawall and adjusts the shoreline 


position as necessary. 


INTEGER YSBEG, YSEND 
REAL Q(41),Y(40),YS(40) , YCOLD (40) 
REAL £(40),F(40),P(41),R(41) 


I=YSBEG 
IF(@(I).GT.0O) TH 


EN 


(Q@) due to limited sand 


Implicit calculation scheme. 


@ positive: Calc of shoreline with correction of @ and Y 


as necessary. 


IF(@(I+1).GE. 


0) THE 


N 


Q(I+1)=P(1+1)#@(1) +R(I+1) 
YC=2¥B¥(Q(1)-Q(I+1))+YCOLD(1) 


IF(YC.LT.Y 


S(1I)) 


THEN 


DQ=(YS(1I)-YCOLD(I))/(2*B) 
Q@(I+1)=Q(I)-Da 


ENDIF 


Y (1) =B¥(Q(I)-Q(I+1))+YCOLD(T) 


IT=I+1 
IF(I.E@.YS 
GOTG 10 
ENDIF 
K=I 
I=I+1 
IF(I.EQ. YSEND 


END+1) 


+1) TH 


GOTO 100 


EN 


Y(I-1)=B¥(@(I-1)-@(I))+YCOLD(I-1) 


GOTO 100 
ENDIF 
IF (1.E@. YSEND 
T=I+1 
GOTO 30 
ENDIF 
ELSE 
K=YSBEG 


) THEN 


Q@ negative: Search for minus piont. 


ENDIF 
IF (Q@(I+1).LT.0O) 
I=I+t1 
IF (1.E@. YSEND 
IF(Q@(I+1). 


THEN 


) THEN 
LE.O) 


THEN 


If absent, calc Y 


Correct @ as necessary. 


YC=2¥B¥(Q(I)-Q@(I+1))+YCOLD(1) 
IF(YC.LT.YS(1I)) THEN 
D@=(YS(1I)-YCOLD(1))/(2%B) 
Q(1)=Q(I+1)+DA 


ENDIF 


Y(1)=B¥(@(I)-Q(I+1))+YCOLD(1) 


GOTO 30 
ENDIF 
ENDIF 
GOTO 20 
ENDIF 


Minus point: Corr of @ out of the element if shoreline moves 


behind seawall. 


YC=2%¥B¥(Q(1)-Q(1+1))+YCOLD(I) 


TFCYE.ET. VS (1)) 


THEN 


AQ 


730 DQ=(YS(1I)-YCOLD(1I))/ (2*B) 


760 QDIFF=Q(1)-@(I+1) 

770 Q(1)=Q(1) ¥DQ/QDIFF 

780 Q(I+1)=Q(1+1)¥D@/QDIFF 
790 ENDIF 

800 Y (I) =BR(Q(1I)-@(I+1))+YCOLD(1) 
810 CX 


820 C¥ Calc of Y starting from element to the left of minus 
830 C¥ point or boundary. @ is negative. 


840 CX 

850 30 DO 40 J=I-1,K,-1 

860 Q(J)=E (J) #Q(J+1) +F (I) 

870 YC=2EBK(Q(J)-Q(J+1))+YCOLD(J) 
880 IF(YC.LT.YS(J).AND.J.GE.YSBEG) THEN 
890 DQ=(YS(J) -YCOLD(J))/(2*B) 
900 Q(J)=Q(J+1)+DAQ 

910 ENDIF 

920 Y(J)=BH(@(J)-Q(J+1))+YCOLD(J) 
930 40 CONTINUE 

940 T=I+1 

930 IF(I.GE.YSEND+1) GOTO 100 

960 CX 


970 C¥ Cele of Y starting from element to the right of minus 
980 C¥ point or boundary. @ is positive. 


990 Ck 

1000 GoTo 10 

1010 100 CONTINUE 

1020 DG 110 I=YSBEG-1,1,-1 

1030 Q(I)=E(1) #@(I+1) +F (1) 

1040 Y(1)=BE(Q(1)-@(1I+1))+YCOLD(T) 
1050 110 CONTINUE 

1060 IF (YSEND.NE.N) THEN 

1070 DO 120 I=YSEND+1,N 

1080 Q(I+1)=P(I+1) #@(1) +R(1I+1) 
1090 Y(1I)=BR(Q(1I)-Q(I+1))+YCOLD(TI) 
1100 120 CONTINUE 

1110 ENDIF 

1120 RETURN 

1130 END 
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SUBROUTINE INDATA(IT,1IT1,1T2,H,Z,N,DT) 
CX SPECIFIES WAVE HEIGHTS AND ANGLES AT SPECIFIED TIME STEPS 


DIMENSION H(40) ,2(40) ,A(40) 
Cc ; 
CXIDUM=0: 
C¥Detached breakwater version. The program gives representative 


C¥wave data (H,Z) simulating effect of shore parallell detached 
C¥breakwater 16*DX offshore and running from I=12 to I=28. 
C¥Inital beach is straight line. 
Cc¥ 
C¥IDUM=1: 
C¥Represents an initially circular beach with no offshore structures. 
CX 
IDUM=1 
DTR=3. 141593/180. 
IF(IT.EQ@.IT1) GOTO 20 
IF(IT.E@.IT2) GOTO 30 
C¥ 
CX¥WAVE ANGLES 
CHXKKXKKHKHKHKKHE 
C¥Case 1: Unaffected breaking angle = 0 deg. 
IF(IDUM.EQ@.1) GOTO 14 
DO 10 I=1,12 
Z(1)=0. 
10 CONTINUE 
DG 11 1=13,19 
Z(I)=Z(I-1)+10./7. 
iba CONTINUE 
Z(20)=-10. 
DO 12 I=21,28 
Z(1)=Z(1I-1)+10./8. 


12 CONTINUE 
DO 13 I=29,N 
Z(1I)=0. 
13 CONTINUE 
GOTO 50 
C¥Alternative case 1: Unaffected angle = -20 deg. 
14 CONTINUE 


DO 15 I=1,N 
Z(1)=-20.¥FLOAT (40-1) /40. 
15 CONTINUE 


GOTO 50 
C¥ 
C¥Case 2: Unaffected breaking angle = -10 deg. 
20 CONTINUE 
IF(IDUM.E@.1) GOTO 27 
DO 21 I=1,8 
Z(1I)=-10. 
21 CONTINUE 


DO 22 1=9,14 
Z(I)=Z(I-1)+10./6. 
22 CONTINUE 
DO 23 I=15,18 
Z(1)=Z(I-1)+2.5 
23 CONTINUE 
DO 24 1=19,20 
Z(I)=Z(I-1)-5. 
24 CONT INUE 
DO 25 1=21,23 
Z(1)=Z(I-1)-2.5 


25 CONTINUE 
DO 26 1=24,40 
Z(1I)=-10. 

26 CONTINUE 
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GOTO 60 


C#Alternative case 2: 


27 CONTINUE 
DO 28 I=1,N 
A(I)=-Z(N+1-I)/ 
28 CONTINUE 
DO 29 I=1,N 
Z(1I)=A(T) 
29 CONTINUE 
GOTO 60 
Cc 
C¥Case 3: Unaffected breaking angle 
30 CONTINUE 
IF(IDUM.EQ@.1) G 
DO 31 I=1,18 
Z(1I)=15. 
31 CONTINUE 
DO 32 I=19,28 
Z(1I)=Z2(1-1)-1.5 
32 CONTINUE 
DO 33 I=29,33 
Z(1)=Z(1-1)+2.5 
33 CONTINUE 
DO 34 1=34,40 
Z(1I)=15. 
34 CONTINUE 
GOTO 70 
C#Alternative case 3: 
35 CONTINUE 
DO 36 I=1,N 
Z(I)=0. 
36 CONTINUE 
GOTQ 70 
C%¥ 
C¥WAVE HEIGHTS 
CHEKHKHHHHKKHKHE 
C¥Case 1 
50 CONTINUE 
IF (IDUM.EQ.1) 
DO 51 I=1,10 
H(I)=1.50 
si CONTINUE 
DO 52 I=11,20 
H(I)=H(I-1)-0. 
52 CONTINUE 
DO 53 1=21,33 
H(I)=H(I-1) +1. 
53 CONTINUE 
DO 54 1=34,40 
H(I)=1.85 
54 CONTINUE 
GOTO 100 
C¥Alternative case 1 
55 CONTINUE 
DO 56 I=1,N 
C% H(1I)=2.5 
56 CONTINUE 
GOTO 100 
Cc 
C¥Case 2 
60 CONT INUE 
IF (IDUM.EQ.1) 
DO 61 I=1,2 
H(I)=1.85 
61 CONTINUE 


H(I)=3.0-FLOAT(I-1)%.5/40. 


Unaffected angle = 


DTR 


oTo 35 


GOTO 55 


1 


35/14. 


GOTO 65 


20 deg. 


= 15 deg. 


Unaffected angle 


Al2 


O deg. 


1400 DO 62 1=3,16 


1410 H(I)=H(I-1)-1.35/14. 
1420 62 CONTINUE 

1430 DO 63 I=17,29 

1440 H(I) =H(I-1)+1.35/14. 
1450 63 CONTINUE 

1460 DO 64 I=30,40 

1470 H(I)=1.85 

1480 64 CONTINUE 

1490 GOTO 100 

1500 C¥Alternative case 2 

1510 65 CONTINUE 

1520 DO 66 I=1,N 

1530 A(I)=H(N+1-I) 

1540 66 CONTINUE 

1550 DO 67 I=1,N 

1560 H(I)=A(T) 

1570 67 CONTINUE 

1580 GOTO 100 

1590 C% 


1600 Cx¥Case 3 
1610 70 CONTINUE 


1620 IF(IDUM.EQ.1) GOTO 75 
1630 DO 71 I=1,12 

1640 H(I)=1.0 

1650 71 CONTINUE 

1660 DO 72 I=13,24 

1670 H(I)=H(I-1)-0.5/12. 

1680 72 CONTINUE 

1690 DO 73 1=25,37 

1700 H(I)=H(I-1)+1.35/14. 

1710 73 CONTINUE 

1720 DO 74 1=38,N 

1730 H(I)=1.85 

1740 74 CONTINUE 

1750 GOTO 100 

1760 C¥Alternative case 3 

1770 75 CONTINUE 

1780 DO 76 I=1,20 

1790 A(I)=H( 21) 

1800 76 CONTINUE 

1810 DO 77 I=21,N 

1820 A(I)=A(N-I+1) 

1830 77 CONTINUE 

1840 DO 78 I=1,N 

1850 H(I)=A(T) 

1860 78 CONT INUE 

1870 C% 

1880 CXOQUTPUT 

1890 CX 

1900 100 CONT INUE 

1910 IHOURS=(IT-1)*INT(DT) 
1920 WRITE (¥, 600) 

1930 IF(IT.EQ.1) WRITE(*¥,%) *INITIAL CONDITIONS’ 
1940 IF(IT.EQ.IT1.OR.1IT.EQ.1IT2) THEN 
1950 WRITE (¥,%¥) *CONDITIONS AFTER’, IHOURS, ’HOURS’ 
1960 ENDIF 

1970 WRITE (*,%¥) *WAVE HEIGHTS’ 
1980 WRITE(¥,602) (H(I),1I=1,N) 
1990 WRITE (¥, 600) 

2000 WRITE(¥,%) ’WAVE ANGLES’ 
2010 WRITE(¥,602) (Z(1),1=1,N) 
2020 WRITE (¥, 600) 

2030 DO 200 I=1,N 

2040 Z(I)=Z(1)¥DTR 


2050 200 CONTINUE 

2060 600 FORMAT (//) 

2070 602 FORMAT ((1X,10F8. 2) ) 
2080 RETURN 

2090 END 

END OF FILE 
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APPENDIX B: NOTATION 


= o 00a 


aw 


Volume of solids/total volume 

Subscript denoting breaking condition 

Wave group velocity, m/s 

Depth of profile closure, m 

Acceleration resulting from gravity, m/s 
Wave height, m 

Subscript denoting position alongshore 


Dimensionless empirical coefficient in the longshore sediment 
transport rate formula 


Superscript denoting time level 

Total number of calculation cells in the model 

Total volumetric longshore sediment transport rate, m3/s 
Conversion factor from RMS to significant wave height 
Stability parameter, m°/s 

Ratio of density of solids to density of water 

Time, s 

Position alongshore, m 


Position on-offshore; shoreline position, m 


Extremal, internally calculated shoreline position, used in the 


implicit numerical solution scheme, m 

Position of seawall on-offshore, m 

Ratio of wave height to water depth at breaking 
Time step, s 

Space step alongshore, m 

Space increment on-offshore, m 

Angle of breaking waves to the shoreline, deg 
Angle of breaking waves to the x-axis, deg 
Superscript denoting a corrected value 


Superscript denoting a quantity at next time step 
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